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Abstract 

Three new models with V-shaped field potentials U are considered: a 
complex scalar field X in 1+1 dimensions with U (X) = \X\, a real scalar 
field $ in 2+1 dimensions with U(3>) = |<&|, and a real scalar field ip in 1+1 
dimensions with U(tp) = ip 8(<p) where © is the step function. Several ex- 
plicit, self-similar solutions are found. They describe interesting dynamical 
processes, for example, 'freezing' a string in a static configuration. 
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1 Introduction 



Lagrangian for a field, let it be a real scalar field 0, usually contains a kinetic 
part, and a potential U (0) which, in the case of physical models, is assumed to 
be at least bounded from below. In most models, however, U(4>) is a smooth 
function of with isolated absolute minima which are reached at finite values of 
0. Furthermore, if O is one of such minima, the second derivative £/"(0 o ) = 
ml > exists and plays the role of length scale (A ~ 1/wq) in the model. 
For a weak classical field such a model is reduced to a free field model, and when 
quantizing it one can apply the harmonic oscillator paradigm with particle creation 
and annihilation operators. Models with these features are very popular - justly 
so because they have plenty of important applications. 

Said above notwithstanding, one can find interesting models which do not fit 
that description. In particular, there exist physically well-motivated models such 
that U" is infinite at the minimum of U [1J. Specifically, the field potential U is 
V-shaped around the minimum, hence the first derivative has a discontinuity at 
= 4>q. In papers flU |2l |3]| we have investigated certain relatively simple models 
of this kind with a single real scalar field in 1+1 dimensions and, moreover, with 
the field potential invariant under the transformation <f>(x) — > —<p(x). Among 
the most interesting findings was a scaling symmetry of the 'on shell' type. This 
symmetry is universal in the sense that it does not depend on neither the number 
of space dimensions nor the number of fields. 

In the present paper we continue the investigations of field theoretic mod- 
els of that kind. Our goal here is to find self-similar solutions in models which 
are less restricted than the ones considered in the previous papers: we allow for 
more fields, or more space dimensions, or smaller symmetry. We consider three 
new models. In the first one we have a complex scalar field in 1+1 dimensions 
(Section 2), the second involves a real scalar field in 2+1 dimensions (Section 3), 
and in the third model the field potential is not invariant under the transformation 
(f)(x) — > —<j>(x) and has degenerate minimum extending from = to — oo (Sec- 
tion 4). These models are interesting already on purely theoretical ground as ex- 
amples of V-shaped field potentials and this is our main motivation. Nevertheless, 
these models can have applications. For example, the first two can be regarded as 
models of pinning of a string or a membrane, respectively, and the third model de- 
scribes the process of depinning of a planar string. It should be stressed that total 
energy is conserved in these models - the pinning occurs because of the dynamics 
and not because of a dissipative loss of energy. 

In the case of V-shaped potentials the pertinent Euler-Lagrange equations are 
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nonlinear in a rather unusual way - as a rule they contain a discontinuous func- 
tion of the field. Restricting the considerations to the sector of self-similar fields 
is a natural simplifying step which reduces by one the number of independent 
variables. Theory and applications of self- similar solutions of nonlinear evolution 
equations is a well-established, important branch of theory of nonlinear systems, 
see, e.g., flH|5]|. It turns out that in the sector of self- similar fields there is quite 
interesting dynamics which can be seen, e.g., from the explicit solutions which 
are presented in subsequent Sections. 



2 Complex scalar field in 1+1 dimensions 

The string in three dimensional space is attracted to a straight line, which we 
call the ^-axis. The two directions perpendicular to it are denoted as the Xi, X 2 
directions. For notational simplicity, we use dimensionless coordinates z, X\ , X 2 - 
the physical coordinates are obtained from them by multiplying by a unit of length. 
Position of the string at a fixed time t is given by two functions X t (z, t),X 2 (z, t) 
of z E (— oo, ooV, which give Cartesian coordinates of points of the string in the 
(Xi,X 2 ) planes [J . We will consider only the cases of a stretched string, so that 
Xl, X 2 are single- valued functions of z. 

We are interested in the dynamics of this system in the special case when the 
attractive force has the potential 

U(X 1 ,X 2 ) = \X\, (1) 

where X = X\ + iX 2 is the complex number representation of points from the 
planes perpendicular to the z-axis. Then, the evolution equation for the string, in 
the approximation discussed in the previous Section, has the form 

(d? - dl)X k = - j^L, (2 ) 

where k = 1,2. Thus, the attractive force has constant modulus equal to +1. 
Equation (2) can be written as 

{%-%)X = F{X), (3) 

where 

nx) = ~. (4) 



Also t is a dimensionless variable proportional to the physical time. 
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Here we have assumed that 1^0. However, one should take into account the 
fact that X = is physically acceptable configuration of the string - the string just 
rests on the z-axis. In order to formally include X = to the set of solutions of 
the evolution equation we assume that F(X) = if X = 0. Of course F(X) is 
discontinuous at X = 0. 

It should be noted that we are looking for so called weak solutions of Eq. (2), 
EQ. When such a solution is inserted in Eq. (2) one should obtain an identity 
only when the both sides are integrated with arbitrary test function of the variables 
t, z, while in the case of strong solutions the identity is obtained immediately after 
the substitution. Actually, the weak solutions are the right ones in the context of 
Euler-Lagrange equations because the stationary action principle has precisely the 
weak form 

/ dtdx 5Xi-^r = 0, 
J oXi 

where S is the action functional and 5Xi arbitrary test functions. 

The potential U is V-shaped: plot of U (Xi,X 2 ) has the form of a symmetric 
cone with the tip at the point X\ = X 2 = 0. Equation (3) possesses the scale 
invariance: if X is a solution of it, then 

X x (z,t) = X 2 X(^ t -) (5) 

is a solution too for any A > OS By definition, the self-similar solutions are 
invariant with respect to these transformations. 
It is convenient to use a polar Ansatz for X : 

X = p exp(zx), 

where p and x are real. Then 

— = sign(p)exp(ix), 

where the signum function takes values ±1, or if p = 0. Eq. (3) is equivalent to 
the following two equations 

p (d 2 t - df) X + 2 (d tX d t p - d zX d zP ) = 0, (6) 

transformations with A < can be obtained by combining transformations (5) with the 
space-time reflection z — > —z, t — > —t. 
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(d 2 t -dfjp- (dtx dtx - d z x d z x) p = -sign(p). (7) 

In the case of constant phase x this set of equations reduces to the signum-Gordon 
equation considered in J3]|. For this reason, we focus here on self-similar solutions 
of Eqs. (6, 7) with non constant x- 

Self-similar Ansatz for X has the form 

X = H(y), p = z 2 S(y),y= t -. (8) 

Then Eq. (6, 7) are reduced to the set of nonlinear ordinary differential equations 
for S and H: 

(1-y 2 ) (SH" + 2H'S') + 2ySH' = 0, (9) 

(1 - y 2 )(S" - H' 2 S) + 2yS' - 2S = sign(S), (10) 

where ' stands for the derivative d/dy. Notice that the presence of the factor 1 — y 2 
implies that values of H' and S' may have finite jumps at y = ±1. One way to 
see it is as follows: the finite jump of the first derivatives would lead to Dirac's 
delta terms ~ 5(y ± 1) in the second derivatives, but such terms are not harmful 
because (1 — y 2 ) 5(y ± 1) = 0. 

Finding solutions with non constant x is greatly simplified by the observation 
that Eq. (9) can be written in the following form 




Hence, 

S 2 H' = c (l-y 2 ), (11) 
with c being a real constant. Therefore, we can eliminate H' from Eq. (10), 

(1 - y 2 )S" - Cl{l ~f ]i + 2yS' - 2S = -sign(S). (12) 

Equation (12) has the following solution, valid in the region \y\ < 1, 

S 1 = a(l-y 2 ), (13) 
where a is a real constant related to cq\ 

4a 4 + c 2 = \a\ 3 . (14) 
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Inserting formula (13) in Eq. (12) and integrating the resulting equation for H we 
find that 

H(y) = ^\n\±y+h , (15) 
2 1-y 

where h is a constant and d = c /a 2 . Formula (14) implies that 

1 



\a\ = 



4 + 



2 ' 



Si = ±r4- (16) 



Therefore, 

i-y 2 

Let us stress that solution (13) is valid only in the interval \y\ < 1. If we take 
\y\ > 1, then instead of (14) we obtain the condition 

4a 4 + Cq = -|a| 3 , 

which implies that a = = cq. Therefore, the Ansatz (13) considered in the 
region \y\ > 1 yields only the trivial solution S = 0. The partial solutions So, Si 
match each other at y = ±1. Together they form the continuous solution on the 
whole y-axis. The first derivative with respect to y has a finite jump at y = ±1, 
but this is allowed for by Eq. (10). 

Inserting z 2 S\ and H in the polar Ansatz we obtain the final form of the self- 
similar solution for t > (at a given time t > the solution Si covers the z > t 
and x < — t parts of the z axis, while the interval \z\ < t is covered by the trivial 
solution S = 0): 



V (^ 2 



-t 2 ) exp (if ln|±| + i/io) if \z\>t, 



x( Z ,t) = v ~ 1 v 2 >-* 1 ^ ' (i7) 

if \z\ < t. 

The constant phase ho is not interesting - it just corresponds to global rotations 
around the z-axis. The minus sign present in formula (16) has been included in 
this phase. It is clear that the distance between the string and the 2-axis, equal to 
\X\, quadratically grows from at the points z = ±t to oo when z — > ±oo. Much 
more interesting is the behaviour of the phase of X at a fixed time t. It describes 
the winding of the string around the 2-axis. At z — > oo the phase \ is equal to h . 
When z decreases, the phase increases by 2n with each step from z to z — Az, 
where 

z 2 - t 2 

Az = 



z + coth(27r/<i )t 
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(we have assumed here that do > 0). Thus, the string winds around the 2-axis 
infinitely many times as z — > t. Similar behaviour of the phase is found when 
z — > —t, the only difference is that the string winds in the opposite direction. 




Such unusual behaviour of the string might create a suspicion that it is un- 
physical for the following reason: that there is an infinite amount of energy accu- 
mulated in the finite region around the point z = t (and, symmetrically, around 
z = —t). It turns out that this is not the case. The total energy E T of the piece 
of the string between the points with the coordinates z = t and z = t + T, where 
T > 0, is given by the formula 

E T = dz {^dtXAXi + ^XidzXi + . 
Elementary calculation shows that 

E ^\ih (T+t)(T+2t) - 
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Furthermore, the length of that piece of the string is also finite: 



t+T . rz z\ tT 



L T (t) = J t dz ^\d z X\ 2 + l 



4 + cZ§ 



The last formula is valid for t > T and t > y 4 + d^. 

Let us notice that if we put Co = 0, then do = 0, and the solution has the simple 
form 

j(z 2 — t 2 ) exp(ih ) for \z\ > t 
for \z\leqt. 



X 



This particular solution coincides with one of self-similar solutions of signum- 
Gordon equation considered in Q. In this case the string lies in one plane con- 
taining the z-axis. Solutions with c ^ are more general - the string winds 
around the 2-axis. In the case of cq = Eq. (12) has the same form as in the 
signum-Gordon model 0. Thus, the 'centrifugal' term ~ Cq in Eq. (12) is the 
only effect of presence of the second scalar field X 2 , as far as the self- similar 
solutions are concerned. 

Solution (17) can be regarded as the solution of an initial value problem for 
the string being pinned to the z-axis. The initial data are given at the time t = 
which corresponds to y = 0. Formula (17) implies that 

X{z,t = 0) = j^eMih o )z 2 , (d t X)(z,t = 0) = ^-^ex P (ih )z. (18) 

Thus, at the initial moment all points of the string lie in one plane (the one 
with azimuthal angle h ), and the initial velocities (d t X)(z, t = 0) of the points 
of the string are perpendicular to their initial position vectors X(z,t = 0) = 
(X 1 ,X 2 )(z,t = 0). 

Let us end this Section with the observation that one can also obtain asymmet- 
ric solutions, i.e., such that X = for z < t or for z > —t. The corresponding 
initial data are obtained by multiplying both formulas (18) by Q(z) or Q(—z). The 
point is that for t > y = corresponds to both z = oo and z = — oo. Therefore, 
Eq. (12) should actually be considered separately in the regions y > and y < 0, 
and there is no justification for a requirement of continuity of S at y = 0!. For 
this reason we may take the trivial solution So in the whole region y < and Si 
in the interval [0,1], or S for all y > and Si for y G [—1,0]. These choices give 
the asymmetric solutions. On the other hand, notice that y = oo and y = — oo 
correspond to z — 0. Therefore, the behaviour of S(y) at y — » oo is correlated 
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with the behaviour at y — > — oo because the function X is continuous at z = 0. 
However, this is not such a severe condition as it might seem because formula (8) 
for the p function contains the factor z 2 which vanishes for z — - it relates only 
the terms ~ y 2 in the function S(y) in the regions y > 1 and y < — 1. 

3 Real scalar field in 2+1 dimensions 

Let us consider a membrane (in three dimensional space) which is attracted to 
a plane. The plane is parametrised by the Cartesian coordinates x, y, and the 
elevation of the membrane at the moment t over the point (x, y) of the plane is 
denoted by <&(x, y, t). We consider the membrane without overhangs, hence <3> is 
a single valued, real function of x, y. The evolution equation has the form 

(d 2 - A)$ = -sign(<f>), (19) 

where A denotes the two-dimensional Laplacian. The corresponding field theo- 
retic potential again has the form (1), that is U (<&) = ]<E»|. The scaling transforma- 
tions which are the symmetry of Eq. (19) have the form 

^(x,y,t) = \ 2 ^,^ t -), 

where A > 0. 

For simplicity, we consider only the axially symmetric configurations, hence 
$ does not depend on the azimuthal angle. In this case the evolution equation has 
the form 

(Of - d 2 r - <9 r /r)$ = -signify, (20) 
where r = \Jx 2 + y 2 . The self-similar Ansatz 

$ = r 2 P(w), w = - (21) 
r 

reduces Eq. (20 ) to ordinary differential equation 

(l-w 2 )P" + 3wP' -4P = -sign(P), (22) 
where ' denotes the derivative d/dw. The corresponding homogeneous equation 

(1 - w 2 )P" + 3wP' - 4P = 0. (23) 
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has the following linearly independent solutions H 



P x = 1 + 2w 2 , (24) 

p _ \ (1 + 2w 2 ) arccosu; — 3wy/l — w 2 if \w\ < 1, 
2 ~ \ -(1 + 2w 2 ) In |w - Vw 2 - 1| - 3wVw 2 - 1 if H > 1. 

From them we can construct partial solutions of Eq. (22) 

P + = a 1 P 1 + a 2 P 2 + \ if P+>0, P_ = /3iPi + /3 2 P 2 -^ if P_<0, (26) 

where «j, are constants. In general, these solutions are valid only in certain 
finite intervals of the w axis, determined by the inequalities P + > 0, P_ < - for 
this reason we have called them the partial solutions. There also exists the trivial 
solution P = (remember that sign(0) = 0), and the static solutions 

P = ±\. (27) 

By glueing together these solutions one can produce various self-similar so- 
lutions of Eq. (22). We will present here just one class of such solutions, which 
describes how the membrane 'freezes' in the static position. These solutions in- 
volve the static solution and the partial solution P + which are glued together on 
the 'light-cone' r = t, i.e., at the point w = 1, see Fig. 2. In accordance with Eq. 
(23), in this case it is sufficient to demand only continuity of P(w). This condition 
has the form P+(w = 1) = 1/4. Simple calculation shows that it is equivalent to 
the condition a.\ = 0. Therefore, the corresponding solution $(r, t) of Eq. (20) 
has the following form 



— if r < t, 

$(r,t) = i [(r 2 + 2t 2 )arccos^-3t v / ^ 3 ^l +f if r > L ^ 



Here t > 0, and we have chosen that branch of arccosu* for which arccos 1 = 0. 
Moreover, 

a > > ~h (29) 



i P 1 has been obtained by inserting a polynomial Ansatz in Eq. (22), and P 2 from Pj by the 
standard method 1 8 1 
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This last condition ensures that P + > on the whole half-line r > t. Thus, with 
time larger and larger part of the membrane becomes static. It follows from the 
explicit form (28) of the solution that at the initial time t — 

$(r, t = 0) = - Ua 2 + -) r 2 , d t $\ t=0 = -Aa 2 r. (30) 

If, for example, the initial slope of the parabola is smaller than 1/4 (a 2 < 0), 
the initial velocity is positive, i.e., the parabola moves upwards towards the static 
shape. Notice, however, that this is done in a manner consistent with the causal 
dynamics implied by wave equation (20) - the freezing occurs at the front r = t 
which travels with the finite velocity equal to +1. 




i i i — i — i — i — i — i — i — i — i — i 1 i — i — i — i — i — i — i — i — i — i — i — i — i 



r = t 

Figure 2: Snapshot of self-similar solution (28). 



One may expect that there are self-similar solutions of Eq. (19) which depend 
on the azimuthal angle 9 in the (x, y) plane. The angle 9 is invariant under the 
rescaling (x,y) — > (x/X,y/X). The Ansatz for self-similar solutions could be 
taken in the form $ = r 2 P(w, 9). We have not investigated such solutions of Eq. 
(19). 
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4 6-Gordon model in 1+1 dimensions 



As pointed out in [1], the scaling symmetry (4) is shared by all (l+l)-dimensional 
models with V-shaped field potential 

U((p) =c 1 (p Q(ip) -c 2 if 9{-<p), 

where ci, c 2 are arbitrary real constants, 6 denotes the step function^, and ip(z, t) 
is a real scalar field. The corresponding wave equations have the form 

(d 2 - d 2 ) V = -cMf) + (31) 

We are interested in the models such that for ip = the potential energy has 
minimum at ip = . For this, we have to assume that c\ > 0, c 2 > 0. In the 
particular case c\ = c 2 = 1 the potential U is symmetric with respect to the re- 
flection ip — > —ip, and Eq. (31) becomes the signum-Gordon equation considered 
in 0. Now we would like to investigate the self- similar solutions in the case 
Ci = 1, c 2 = 0, when Eq. (31) has the form 

{d 2 t -d 2 z )^ = -Q{ V ). (32) 

For the obvious reason we will call this equation the ©-Gordon equation. 

Physical context for considering equation of the form (32) is the depinning 
phenomenon. This equation can be regarded as describing the dynamics of a 
string in a plane which would be permanently pinned to the z line were it not for 
a constant bias force, which exactly compensates the pinning force on one side of 
the z line (ip is just the deviation of the string from this line). Hence, the 6-Gordon 
equation describes the dynamics of the string exactly at the depinning transition. 

The Anstaz for self- similar solutions has the form 

(p(z,t) = z 2 R(y), 

where y = t/z. Inserting it in Eq. (32) we obtain the following equation 

(1 - y 2 )R" + 2yR' - 2R = -&(R)Q(z 2 ) 

where ' denotes the derivative d/dy. This equation is to be satisfied in the weak 
sense, hence we may drop the factor Q(z 2 ) which differs from 1 only at the single 

4 Q((p) — 1 when ip > and Q((p) — when ip < 0. 
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point 2 = 0. Therefore, R(y) obeys (again in the weak sense) the following 
equation 

(1 -y 2 )R" + 2yR'-2R = -Q(R). (33) 

Equation (33) has the trivial solutions R = 0, R = 1/2, which correspond, 
respectively, to ip = and to the static solution 

Vs = z 2 /2. (34) 

Furthermore, there are the partial solutions 

R+ = a 1 + a 2 y+(a 1 -^)y 2 if R + > 0, R- = h(l + y 2 ) + b 2 y if i?_ < 0, 

(35) 

where <2j, hi are constants. As in the previous Sections, the partial solutions are 
valid only in certain intervals of the y axis, determined by the inequalities R + > 

;J R_ < 0. 

Glueing together the partial solutions one can obtain whole variety of self- 
similar solutions of Eq. (32) [9]. We present here certain interesting examples. 
1. Freezing in the static configuration. 
This solution is depicted in Fig. 3. 
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It has the following explicit form 



(f(z,t) 







±z 2 

2^ 



2 {V1V2 —) (Vxt ~z){z- V 2 t) 



z<0, 
ze[0,t), 
z e 

Z > V\t, 



(36) 



where 



2|5 | 



> 1, v 2 



Vi 



2v 1 



< 1, 



and S < 0. The parameters S^, S are not independent - they have to obey the 
equation 

S'J + (2S - 1)S' 2 - 8S 3 + AS 2 -So- AS 2 S' = 0, (37) 

which follows from the condition of continuity of the derivative d z ip at the point 
z = Vit. It turns out that cubic equation (37) has three real solutions when S < 0. 
The relevant root S' obeys the inequality S' > 1/2 — 2S . 
Initial configuration for solution (36) has the form 

s z 2 e{z) 



d t <p = S'qzQ(z). 



2. Depinning of the string from the z-axis. 

The function <p has the shape presented in Fig. 4. 




Figure 4: Snapshot of self-similar solution (38). 
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In this case 



f 



(p(z,t) 



where Sn > 0, and 



-a(z + t) 2 



[ S z 2 + S' tz + (S - \)t 2 



Vi 



S' + y/^ 2 - AS 2 + 2S 
2S n 



a = 1(250^ + ^)^^-1 <1, 
4 v i + 1 



2 < ~t, 

z e [-t,t], 

z E [t,Vit], 

Z > fit, 



>1, 



(38) 



(«i - I) 2 

Solution (38) is valid provided that > 1. This condition is satisfied if 



< So, S' < 1/2 - 2S . 



(39) 



5 Remarks 

1. We have presented examples of self- similar solutions which describe rather 
interesting dynamical processes like the freezing in the static configurations or 
the depinning. However, it is obvious that the three models considered here have 
other self-similar solutions as well. In the case of 6-Gordon model one can pro- 
vide the complete list of such solutions |[9l analogous to the one given in for 
signum-Gordon model. In the other two models such a complete list is probably 
out of our reach because of the S~ 3 term in Eq. (12), or possible azimuthal angle 
dependence in the case of membrane model. 

2. The kinetic part of evolution equation (3) contains d'Alembert operator d 2 — d 2 , 
and this is standard for bosonic field. On the other hand, when we apply this 
equation in order to describe evolution of the string we automatically make certain 
assumptions about the kind of string we consider. The precise statement is that it 
is the string which has evolution equation of the form (3) when its world-sheet is 
parametrized by the coordinates t, z. Of course, we would like to see a connection 
with Nambu-Goto string, which has evolution equation of the form 

9 ' ^9'^) = 0, (40) 



du a \ du b 
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where u° = t,u 3 = z,(x^) = (t, Xx{z,t), X 2 (z,t), z), (g ab ) is the inverse to 
(gab), and 

dx^ dx v 

Let us take a slightly more general than Minkowskian space-time metric r/^: 

(^)=diag(l,-/2,-/ 2 ,-l), (41) 

where Zq is a constant (equal to 1 in the case of Minkowski space-time). Simple 
calculations show that Nambu-Goto equation (40) is reduced to the equation 

dfX { - d 2 z Xi = 

when 

lodtXi^l, kdzXi^l. (42) 

Interesting possibility to satisfy the conditions (42) is to take a very small l . This 
would correspond to the Nambu-Goto string (or a linear defect like a vortex) in an 
anisotropic medium which would effectively provide metric (41). The term F(X) 
in Eq. (3) represents the pinning force with which the z line attracts the string. 

Analogous remarks can be made about the membrane discussed in Section 3 
and about the planar string of Section 4. 

3. The V-shaped field potentials U(<f>) we consider should not be identified as 
potentials which just contain the modulus of the pertinent field. Examples of such 
potentials can be found in |[T0~l[TTll . The fundamental difference is that in our mod- 
els the second derivative U" does not exist right at the minimum of the potential, 
hence there is no preferred finite length (or mass) scale. Potentials considered in 
IfTOl [TTi have the A shape. In these cases U" does not exist at a local maximum 
while at minima it exists - that has much smaller impact on properties of the fields. 
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